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1. Introduction 

Consider a system of N charges in a spherical container of radius R. An external 
potential centered at the origin exerts an attractive central force. The Hamiltonian is 




where r^^ and are the position and velocity of charge a. The repulsive interaction 
potential between particles a and 7 is denoted by V^Va-y) where r^-y = Ir^ — r^|. The 
external "confinement" potential Vc (r^) is the same for all particles, and Vw {ra) is the 
wall potential that is zero inside the container and infinite otherwise. Charged systems 
are of direct relevance e.g. to dusty plasmas and ions in traps. For neutral systems 
Vb (ra) is the interaction of each particle with a uniform neutralizing background, 
corresponding to an OCP or jellium with an attractive trap at the origin. 

An example that has been studied recently is jellium with a point positive ion of 
charge number Z at the origin [H [2] • At equilibrium the electron density is enhanced 
(partial confinement) near the origin and approaches a uniform limit far from the ion, for 
sufficiently large R, due to charge neutrality. An opposite extreme is a system of charges 
in a strong trap such that they are localized in a finite domain away from the wall and 
their density vanishes outside this domain (complete confinement). There is a continuous 
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mapping between these two limiting cases, controlled by the relative strengths of the 
repulsion between particles and the strength and form of the confining potential. The 
discussion of these cases here will be limited to classical statistical mechanics. 



\nn{ri) + pVc(ri) - / rfrg (n (rs) - ri;,) c (rai) 



0. (3) 



2. Equilibrium density 

In the absence of a confining potential the equilibrium density may be uniform (charge 
neutral) or nonuniform (charged). In general the confining potential will induce a 
nonuniform equilibrium density in any case. Formally it is determined from the exact 
Yvon-Born- Green equation 

dn(ri) ^ , . dVc(ri) ^ f , , , . , , ^ dV (r2i) , . 

where ?T,(ri,r2) is the joint density for two charges and n\, is the density of the 
uniform neutralizing background. A reasonable approximation is the hypernetted chain 
approximation (HNC) [3lHj, 
d_ 

dri 

The solutions to this equation are quite rich, reflecting the competition between the 
attraction of Vc and the renormalized repulsion represented by c(r2i). They have been 
described in some detail for the case of the positive ion in jellium [1]. Here we report on 
some preliminary results for the qualitatively different case of weakly coupled charges 
g in a harmonic trap Vf. (r) = kr"^ /2 with U}, = 0. There is then a competition between 
Coulomb repulsion and attraction by the central trap, with charge density enhanced 
near the wall or near the center, respectively. These effects are exactly balanced when 
uim = ujc, leading to a uniform density at all temperatures. Here tUm = ujp/ y/S is the Mie 
plasma frequency = Airnq'^/m) and Uc = \fkfm is the frequency associated with the 
harmonic trap. More generally, when lOc > Wm the charges are increasingly drawn away 
from the wall as the temperature is lowered, approaching at T = a uniform density 
for r < Rq = {ijjm/ujc)^^^ R and zero density for _Ro < < -R 0. In the opposite case of 
ujc < uJm the Coulomb repulsion dominates and the particle density is enhanced at the 
wall. As the temperature is lowered the density becomes sharp at the walls - a "Coulomb 
explosion" that is restrained by the external walls. This behavior is illustrated in Fig. 
la) for uoc/ujm = 2 and in Fig. lb) for Uc/wm = 1/2, at T = 1, 0.1, and 0.01. 



3. Linear Response 

Consider a perturbation of the nonuniform equilibrium state by an external potential 
of the form Uextit) = Y.a=i^ext {ra,t) = f drVext{r,t)n{r), where n(r) is the phase 
function representing the particle density. Then, to linear order in Ve^t {f^, t) the response 
of the average particle density to this perturbation is [6] 

S {n{r,t)) = j^dt' j dr'x (r, r', t - t') Kx* (r', t') , (4) 
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Figure 1. Temperature dependence of densities for trap-dominated (left) and 
Coulomb-dominated (right) conditions. 



with the response function 

X (r, r', t) = -(3dt {n (r) n (r', -t))^ = /5Vr' ■ (n (r, t) j (r')), • (5) 

The brackets ()g denote an equihbrium ensemble average. The second equahty follows 
from the conservation law dtfi (r, t) + Vr ■ j (r, t) = 0, and the stationarity of the 
equilibrium ensemble. The response function x (r, r', t) provides an appropriate object 
for studying the possible modes of excitation that are supported by the system. An 
alternative equivalent instrument is the dielectric function e (r, r'; t) defined in analogy 
with electrodynamics 

e-' (r, r'; t) ^ S (v - r') S (t) + j dr"x (r, r"; t) V(|r" - r'|), (6) 
4. Markov Kinetic Theory 

The evaluation of the response function is a difficult many-body problem. The conditions 
of interest include both strong confinement and possibly strong charge correlations, so 
there is no small parameter available for simplifications. Instead, a non perturbative 
Markov kinetic theory has been discussed recently in this context [2J. It is based on 
approximating the formal generator for dynamics in the single particle phase space by 
its exact form at t = 0. This leads to a mean field theory of the linear Vlasov form, 
but with both the confining potential and the charge - charge potential renormalized by 
the initial equilibrium correlations. The analysis of reference [2] can be extended in a 
straightforward way to the system considered here with the result 

X (r, r', t)^ dt' J dv{e-^ (r, n; t - t') Xo (ri, r', t') . (7) 

where Xo (i""; i"'; ^) is the response function for confined particles without the interparticle 
interactions (V^(|r — r'|) = 0) 

Xo (r, r', t) = -I3n (r) / rfv0 {v) e"^«V • Vr5 (r - r') , (8) 
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(f ) is the Maxwellian, and where the generator Cq for the effective single particle 
dynamics is 



The renormalized confinement potential is determined from the equilibrium density by 
Vc (r) = —f3~^ In n (r) . Also, e (r, r'; t) is 



where the renormalized charge - charge potential, V(r, r'), is defined in terms of the 
direct correlation function by V (r, r') = — /9~^c (r, r'). 

Further progress requires evaluation of xo (r, r', t) whose dynamics is determined 
from the effective single particle charge dynamics in the presence of the confinement 
potential. This can be quite complex over the whole range of weak to strong confinement, 
and has been discussed in some detail for the case of a central positive ion in jellium [2]. 
The autocorrelation function for the total force on the charges was studied as a function 
of the charge number on the central ion. It was found that a simple representation of 
the single particle trajectories in terms of bound and free states provided an accurate 
analytical description from weak to strong ion - electron interaction. 

Preliminary studies of the harmonic trap suggest the possibility of a simplification 
in that well. As noted above, the low temperature limit leads to a uniform 

density inside a sphere of radius Rq. Since the effective potential Vc {r) for the dynamics 
is determined from this potential, xo (r, r', t) becomes the correlation function for free 
particles inside a sphere. An exact evaluation of the force autocorrelation function in 
terms of all response frequencies is possible in this case [7]. 
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